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Background 



• To-date CFD analysis has been performed on, 2D cross- 

sections, 3D extrusions of 2D cross-sections, and 3D ice shapes 


generated by ice accretion codes 



www.nasa.gov 



National Aeronautics and Space Administration 


Motivation 

• Complex 3D ice shape geometry data can now be 
collected 

- In-situ laser scans of ice accretion shapes 

- CAT scans have also been performed 

- Complete ice shape documentation, including surface 
roughness elements 

• How good is good enough? 

- What level of ice shape detail must be simulated by ice 
accretion codes? 

- Detailed analysis of the aerodynamics and heat transfer 
mechanisms at the ice-liquid-air interface can shed light on 
the parameters of importance 
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Ice Accretion Shapes 




• Types of ice accretion 

- Rime 
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Ice Accretion Shapes 
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Workflow 



Ice Shape 
Scanned in IRT 



Ice Shape Processed 
with Geomagic 



Ice Shape Grid 
Built in Pointwise 



Aerodynamic 
Analysis with NCC 



Plotting and Visualization 
with Tecplot 
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Approach (Grid Generation) 



Geomaqic 

Commercial software 
used to create 
watertight surface 
from scanned point 
cloud data 



Lee, S., Broeren, A. P., Addy, H. E., Jr., Sills, R., and Pifer, E. M., “Development of 3-D Ice Accretion Measurement Method,” 
NASA/TM-20 12-2 17702, AIAA Paper-2012-2938, 2012 
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Approach (Grid Generation) 



Pointwise 

Commercial software used to import ice shape geometry data and 
create grid for CFD analysis 

1. Import Geometry 

• Database 

• Surface Grid 
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Approach (Grid Generation) 



Pointwise 

1. Import Geometry 

2. Create Surface Grid - Rime 
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Approach (Grid Generation) 

Pointwise 

1. Import Geometry 

2. Create Surface Grid - Horn 
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Approach (Grid Generation) 

Pointwise 

1. Import Geometry 

2. Create Surface Grid 

3. Create Volume Grid - Rime 
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Approach (Grid Generation) 



Pointwise 

1. Import Geometry 

2. Create Surface Grid 

3. Create Volume Grid - Horn 
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Approach (Grid Generation) 



Pointwise 

1. Import Geometry 

2. Create Surface Grid 

3. Create Volume Grid 

4. Refinement 
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Statistics of Initial Grids 



Ice Shape Chord length Span length Grid Type Volume grid 
Geometry (in) (in) cell count 


Clean - 18 12 Structured 0.5 million 

Rime ED1966 18 6 Unstructured 1.6 million 

Glaze EDI 978 18 6 Unstructured 3.7 million 

Broeren, A.P., Addy, H.E., Lee, S., and Monastero, M.C., “Validation of 3-D Ice Accretion Measurement Methodology for Experimental 
Aerodynamic Simulation,” AIAA 6 th Atmospheric and Space Environments Conference, Atlanta, GA, June 16-20, 2014 




EDI 966 


EDI 978 
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National Combustion Code (NCC) 



• Solver 

- Finite-volume 

- Explicit, four-stage Runge-Kutta integration algorithm 

- RANS, URANS 

• Turbulence 

- k — e model 

- higher order, non-linear method 

- Partially Resolved Numerical Simulation (PRNS) 

• Parallel Computing 

- Parallel Virtual Machine (PVM) 

- Message Passing Interface (MPI) 


Liu, N.-S. and Shih, T.-H., “Turbulent Modeling for Very Large-Eddy Simulation,” AIAA Journal, Vol. 44, No. 4, April 2006 
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Domain Conditions 



• Boundary Conditions 

- Velocity Inlet 

- Pressure Outlet 

- No-slip Airfoil Wall 

- Periodic Side Walls 


Pressure Outlet 



Pressure 

Outlet 


• Freestream Conditions 
M = 0.10, 0.18 
Re = 1 .0x1 0 6 , 1.8x10 6 
P M = 98,595 [Pa] 

T. = 294.3 [K] 
a = 0° to 10° 


Outlet 


Periodic 


Periodic 


Outlet 


Front View 
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Clean Wing (M=0.10@0°) 
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Clean Wing (M=0.10 @ 0°) 




* * * 



* 

#/ r '$ 
■ 

l 

** * * 

i i 

i i 

Computational Results 

* • Experimental Results 




f>. « 

m 



■ * * Zt 

.... 






* 

** 

i 






i 

d 

4 

- 

ji 

1 

► 





1 

i 

i 

i ! 

i i 

' i 



1.5 1 1 1 1 1 

— 0,2 0,0 0.2 0.4 0.6 0.8 1.0 

Normalized Chord, (x/c) 


www.nasa.gov 



Pressure Coefficient, C\ 


National Aeronautics and Space Administration 


Clean Wing (M=0.10@0°) 

Other CFD Solvers 
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Clean Wing (M=0.10@0°) 

Other CFD Solvers 
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Clean Wing C L Curve (M=0.10) 
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ED1966 Wing (M=0.10 @ 0°) 

Rime Shape 
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ED1966 Wing (M=0.10 @ 0°) 

Rime Shape 



- 1.0 



— 
* __ 

i 

! 1 

i i 

Computational Results 

• * Experimental Results 



ft 

* */ 

* : 

* * ft ^ , 

-r" t 

-• *:■*, . .* * ■ 

- 7 r,\ ~ *• 

: - f r ; .^ ■ r-r . 

-L 



tl _ • * < 

.* * ^ 

‘" 4, ' ,,r 

L, 

■■ 



h 

£ T 





*1 

'» 

J. 

• l 






L 1 

L 

i j 

1 i 

L i 

i 3 

i 


- 0.5 


QQ 


0.5 


1.0 


^ 0.2 


0.0 


0.2 0.4 0.5 

Normalized Chord, ( x/c ) 


0.8 


1.0 


www.nasa.gov 



National Aeronautics and Space Administration 

EDI 966 Wing Lift Coefficient Results 

Rime Shape 
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Re = 1.0x10 s , M = 0.1 


Re = 1.8x10 6 , M = 0.18 


• Results suggest that viscous effects play a role for the rime ice case, 
consistent with expectations 

• Results from a single instantaneous pressure profile, used in the computation, 
need to be replaced with time averaged and spatially integrated results 
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ED1978 Wing (M=0.18 @ 0°) 

Glaze shape 
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ED1978 Wing (M=0.18 @ 0°) 

Glaze Shape 
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Future Work 

• Detailed examination of solutions 

- Both ice shapes (EDI 966 and EDI 978) 

■ Variations in flow field results across the span 

■ Time averaging of unsteady results 

■ Spatial integration across the span 

■ Grid resolution studies 

■ Turbulence models 

- Glaze ice shape (EDI 978) 

■ Investigate cause of pressure fluctuations near leading edge 

• Parametric study of mesh quality 

- Establish minimum amount of grid points along airfoil surface 

• Perform detailed analysis of ice surface roughness region 

• Develop post-processing modules for NCC to calculate 
standard external aerodynamic parameters 
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Thank You! 


Questions? 


